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The Method of Increments (Mol) has been employed using a multireference approach to calculate the disso¬ 
ciation curve of beryllium ring-shaped clusters Be ra of different sizes. Benchmarks obtained through different 
single and multireference methods including the ab initio Density Matrix Renormalization Group (DMRG) 
were used to verify the validity of the Mol truncation which showed a reliable behavior for the whole disso¬ 
ciation curve. Moreover we investigated the size dependence of the correlation energy at different distances 
in order to extrapolate the values for the periodic chain and to discuss the transition from a metal-like to a 
insulating-like behavior of the wave function through quantum chemical considerations. 


I. INTRODUCTION 


Metal-insulator transitions (MIT)i - — in realistic sys¬ 
tems are often considered in a one-particle picture to 
explain the change of conductivity in a material. In 
situations such as band, Peierls and Anderson insula¬ 
tors, this yields a successful description, but in the case 
of Mott insulators the electron-electron correlation gives 
rise to a contribution which is more important than the 
electron-ion interaction in the localization of the elec¬ 
tronic wave function (WF)&2. In this situations a proper 
description of the correlated WF becomes crucial, which 
cannot be in general achieved through the first princi¬ 
ple methods usually employed for periodic systems such 
as the Density Functional Theory (DFT)£~— . Among 
the different wave function-based local approaches em¬ 
ployed to overcome the dimensionality problem in cor¬ 
related extended and periodic systems, the method of 
increments (Mol)^. — has gained a particular attention 
in the last decade. This approach can be applied to any 
wave function method and as shown in the recent work 
by Voloshina and Paulus^ it successfully retrieves al¬ 
most 100% of the correlation energy of bulk materials 
with large static correlation if applied in a multirefer¬ 
ence (MR) fashion. 

In our previous investigation^ we have exploited the 
quantum chemical version^Sr— of the density matrix 
renormalization group (DMRG) 25 ’ 26 approach to calcu¬ 
late the ground state energy of a model system, i.e. 
beryllium ring-shaped cluster, and explored the use of 
quantum information theory (QIT)2£ — to characterize 
the wave function and thus determine the metal-like and 
insulating-like character of a system in different regions 
of the potential energy surface (PES). In the present pa¬ 
per we will focus again on the same system exploring 
the use of the method of increments (Mol) for closed- 
shell systems obtaining whole ground state dissociation 
curves through a multireference approach that allows to 


describe the crossing region where single reference ap¬ 
proaches such as CCSD(T)^r— fail. DMRG calculations 
as described in our previous work will be used as refer¬ 
ence for testing the Mol approximation to the correlation 
energy. 

Even though a model system is considered, we will un¬ 
derline how the use of standard canonical methods be¬ 
come prohibitive because of the high correlation effects 
involved especially if aiming to the thermodynamic limit. 
As we will show the Mol formalism used in this work in¬ 
volves only localized orbitals (LOs) which allows it to use 
limited active space despite large systems are involved. 
We will exploit this tool to calculate the correlation en¬ 
ergy of Be„ rings as large as n = 90 and extrapolate the 
behavior at the thermodynamic limit. 

This paper is structured as follows: in section EH we de¬ 
scribe the formalism of the method of increments for 
single- and multireference methods as applied in this 
work; in section m the problematics in describing the 
system with canonical wave function methods are un¬ 
derlined and we give the details of our calculations; in 
section m we focus on th results obtained for the Beg 
ring in order to highlight the accuracy and the advan¬ 
tages of the multireference Mol in comparison different 
methods; we also report and compare the results obtained 
for larger rings and the behavior at the thermodynamic 
limit is extrapolated; our conclusions are finally drawn 
in section El 


II. THE METHOD OF INCREMENTS 
A. General Formalism 

The method of increments exploits the short range na¬ 
ture of the electronic correlation. Within this approach 
localized orbitals are used in order to describe the corre¬ 
lation energy as sum of individual contributions coming 
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from the correlation of different parts of the systems, to 
which we will refer as bodies. In its general formalism 
which employs single reference methods such as coupled 
cluster (CC) or perturbation theory approaches, only the 
occupied orbitals are localized and some of them together 
with the virtual canonical orbitals are used to build the 
correlation space. 

In the Mol, one starts from a first approximation of the 
correlation energy, El orr , which is given by the sum of 
all independent contributions each arising from the i th 
body: 


i 

We refer to the individual terms e,; as 1-body increments. 
The contribution expressed in Eq. [T] ranges typically be¬ 
tween 60% and 90% of the correlation energy achievable 
through the chosen correlation method. The remaining 
part of E COII is enclosed in the higher order increments 
which consider the correlation among many bodies. The 
second natural step is then to include one more body and 
define the 2-body increments Aas: 

Aejj = Cij — (€i + 6j ) (2) 

Going forward to higher order increments we get simi¬ 
lar expressions. For instance the 3-body increments are 
calculated using the expression: 


A £ijk — f-ijk — (A Cij + A Cjk + A 6ik) — (ei + + £fe) (3) 

Combining these contributions, one can finally express 
the correlation energy E corl : 

Ecorr - 5 Z €i + Ae b + - ( 4 ) 

i i<j i<j<k 

As pointed out above, the electronic correlation is in gen¬ 
eral short ranged, and the ee-interaction is a two-particle 
interaction which implies that the increments decrease 
with the distance and the order of the bodies, i.e. if the 
following convergence criteria are fulfilled: 


1A Cij 

I > |Ae*ife| for nj > r ik 

(5) 

1A £ij 1 


(6) 


This allows to truncate the expansion of Eq. [H yield¬ 
ing a meaningful amount of the correlation energy. The 
local nature of the Mol and the possibility of truncat¬ 
ing Eq. [4] are crucial points in making the method a 
candidate for the application to extended and periodic 
systems, for which canonical wave function methods are 
generally prohibitive. Nevertheless the Mol cannot be 
universally applied and fails when the required conver¬ 
gences explicated in Eq. [5] and [5] do not occur. Moreover 
since the increments might have alternating sign, the Mol 
is not a variational method. 


B. Multireference Formalism 


As recently shown by Voloshina and Paulusi^, a mul¬ 
tireference Mol approach can successfully be used to cal¬ 
culate the cohesive energy of bulk metals for systems 
with a high static correlation contribution. In compar¬ 
ison with the single reference formalism, a different lo¬ 
calization pattern is required in this approach. Indeed, 
besides the occupied orbitals, also the virtual orbitals im¬ 
portant for the evaluation of the static correlation have to 
be localized. This way one can then calculate incremen¬ 
tal static contributions to E COTT arising from Complete 
Active Space Self Consistent Field (CAS-SCF)^ calcu¬ 
lations performed within the bodies constituted by occu¬ 
pied and virtual LOs. Finally on top of the CAS-SCF 
wave function, a MR calculation is performed for each 
term of the incremental expansion, including the remain¬ 
ing delocalized virtual orbitals. The scheme is sketched in 
Fig. [T] The equations described in the previous section 
remain still valid are used to expand and truncate the 
correlation energy. As described in the following section 
we will employ a minimal basis set and localize all vir¬ 
tual orbitals which will be use to construct the complete 
active space. We will not perform any multireference cal¬ 
culation on top of the CAS-SCF wave functions and we 
will refer to the method as CAS-MoI. 
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FIG. 1. (Color online) Sketch of the partition of the active 
orbitals for the method of increments. In its multireference 
formalism, CAS-SCF calculations are performed using local¬ 
ized orbitals only as shown in the highlighted boxes in order 
to obtain the non-dynamical incremental contributions to the 
correlation energy. On top of the obtained CAS wave func¬ 
tions, multireference calculations can be performed including 
excitation to the delocalized virtuals. 
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III. CALCULATION DETAILS 
A. The System 

As subject of our investigation we have chosen a 
model system characterized by a high static correlation 
contribution, i.e. beryllium ring-shaped clusters, whose 
pseudo-onedimensional structure was chosen to resemble 
the periodic Born-von-Karman boundary conditions and 
to obtain a WF that will inevitably converge with size to¬ 
wards the thermodynamic limit of an infinite linear chain. 
Because of the quasi degeneracy of valence 2s and virtual 
2 p orbitals (more than 93% of the correlation energy of 
Be atom in its 1 S ground state is static), 4n active or¬ 
bitals would be required to obtain a size-consistent CAS- 
SCF reference for MR calculations for Be n . This leads 
to CAS(2n, 4 n) calculations which become of course pro¬ 
hibitive as the number of Be atoms increases. We report 
in Fig. [2] the dissociation curves of Beio calculated using 
a minimal atomic basis set and different active spaces 
within the CAS-SCF method. As one can see, despite 
the use of larger and larger active spaces gives a finer 
and finer description of the potential energy curves, the 
calculated dissociation plateaus he at much higher en¬ 
ergies than the dissociation limit calculated within the 
same method using an active space consisting of the 2s 
and 2 p functions of the free atom. Similar results are 
reported in more details in Table U as described later on. 
This model system is complicated further by the fact that 
the p functions, besides enlarging the active space, play 
a big role at the Hartree-Fock level too. Indeed, while at 
dissociation the HF orbitals will be linear combination 
of pure 2s atomic orbitals, as the interatomic distance 
shortens their p-character increases and, at a certain ge¬ 
ometry, the HOMO switches from a pure s to a pure p 
molecular orbital. In other words, we encounter a cross¬ 
ing between two HF configurations which dominate the 
ground state at different interactomic distances. This 
holds for the finite pseudo-onedimensional clusters in- 
depentely from their size as well for the periodic chain 
given that both the energy and the character of the dis¬ 
crete Hartree-Fock molecular orbitals of Be n rings con¬ 
verge towards the crystal orbitals. In order to supplye a 
graphical representation of these HF wave functions, we 
report in Fig. [3] the band structure calculated in the two 
regimes and the s-character of the Hartree-Fock valence 
orbitals of a periodic beryllium chain close to the min¬ 
imum of the dissociation curve (blue line) and towards 
dissociation (red line). Respectively, we will refer to these 
as configuration (Conf) 1 and 2 and we will indicate the 
two regimes where each of them dominates as metal-like 
and insulator-like regime. Clearly an accurate single ref¬ 
erence method such as CCSD(T) will work perfectly in 
the regimes where one of the configuration is dominant, 
but not around the crossing region. We will show how 
the size problem and the necessity for a MR approach 
will be overcome through the use of the multireference 
method of increments allowing us to describe the PES of 


Be n till dissociation. 

By unitary transformation of the canonical orbitals of 
the two main HF configurations different sets of LOs 
are obtained (see Fig. [3]). Indeed, in the insulator-like 
regime localized orbitals resemble atomic 2s and 2 p or¬ 
bitals, while in the metallic regime a g - and er„-like or¬ 
bitals appear. As we will show, the use of the different 
starting Hartree-Fock configuration has a huge impact on 
the effectiveness of the method, despite in both cases the 
result should converge toward the Full-CI limit. 


B. Computation and basis set 

The MOLPRO quantum chemistry package^ was used 
to perform the different steps of Mol calculations i.e. 
Hartree-Fock calculations, Foster-Boys localization^ and 
CAS-SCF calculations. We employed a minimal basis set 
consisting of Is, 2s and 2 p atomic functions derived by 
the Dunning’s cc-pVDZii. We kept the core Is orbitals 
frozen during the incremental calculations, focusing on 
the correlation of valence orbitals. The Crystal09 42 code 
was employed to calculate the Hartree-Fock wave func¬ 
tion of the periodic chain using the same basis set. 
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FIG. 2. (Color online) Dissociation curves calculated for the 
ground state of Beio at the Complete Active Space Self Con¬ 
sistent Field (CAS-SCF) level of theory, using different active 
spaces. The dashed line corresponds to the dissociation limit 
calculated using a CAS(2,4) for the Be atom. A minimal basis 
set has been used in all cases. 


IV. RESULTS 
A. Be6 ring 

We start our analysis reporting the single reference 
CCSD(T) results obtained for Be6 using both the canon¬ 
ical formalism and the Mol expansion, to illustrate how 
within the same quantum chemical method, this trun¬ 
cation works in retrieving the correlation energy. The 
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FIG. 3. (Color online) Hartree-Foc.k band structure and s- 
character of the valence crystal orbitals of an infinite beryl¬ 
lium chain for the two leading configurations Conf 1 (blue 
line) and Conf 2 (red line). The s-character, defined as the 
normalized sum of the squared coefficients of the s-bases, is 
reported in reciprocal space as well as the valence band for the 
two regimes. Notice that the character of the valence band 
(full lines) at the X-point of the Brilluoin zone switches from 
pure p to pure s as sketched in the insets. The localized or¬ 
bitals obtained by unitary transformation of the valence and 
virtual orbitals are also shown. As one can see, the difference 
in s-character is reflected by the different symmetry of the 
localized orbitals. 


method of increments at the 4-body level retrieves more 
than the 99.99% of the correlation energy achievable with 
the canonical CCSD(T) if the proper HF starting con¬ 
figurations (Conf 1 and 2) are used in the two limiting 
regimes, while around the crossing this percentage drops 
to around 95%. Also in this regime the 7i diagnostic 
calculated within the canonical CCSD(T) presents val¬ 
ues larger than 0.03 highlighting the necessity of using a 
multireference approach for describing the static corre¬ 
lation. However as already stated and shown in Fig. [2j 
canonical CAS-SCF approaches cannot be used to calcu¬ 
late a size consistent wave function unless a CAS(2n,4n) 
is performed, i.e. a Full-CI in our minimal basis set, 
and this is of course not applicable as n increases. Of 
course, a multireference calculation performed on top of 
an accurate CAS or RAS calculation would reach the re¬ 
quired accuracy, but since it is our interest to explore the 
convergence toward the thermodynamic limit, it is clear 


that a local method is preferable. Moreover, this choice 
in general reduce the problem of choosing a proper and 
consistent active space. 

In Fig. 0] we report the dissociation curve of Be6 as 
obtained using the CAS-Mol at the 4-body level using 
the Conf 1 and 2. Because in our case, the incremen¬ 
tal scheme is converging towards the Full-CI solution, it 
would be in principle equivalent to start from one or the 
other configuration. However as it can be seen, truncat¬ 
ing the correlation energy at the 4-body level the two 
results match in a narrow regime only, i.e. around the 
crossing, while for other internuclear distances a non¬ 
monotonic behavior is obtained if not the proper HF 
configuration is used. We want to focus now on the com¬ 
parison with the DMRG data which were obtained in 
our previous investigation^ and we use fixed M = 1024 
block states results as a reference to evaluate the behav¬ 
ior of the CAS-Mol and the other methods discussed. In 
Fig. [5] we report the absolute error with respect to the 
DMRG data of the above described CAS-Mol, canoni¬ 
cal CCSD(T) and the very accurate multiconfigurational 
calculation, RAS (4,24). As it can be seen the CAS-Mol 
gives the best results for almost any interatomic distance, 
and while the other approaches depart too much from 
the DMRG reference around the crossing, the CAS-Mol 
results stay within an error of ± 4 m Eh on the total 
energy. Moreover since the increments are all negative 
and converging (see Fig. [G]) we might argue that the Mol 
is behaving variationally in our calculations. As a conse¬ 
quence, if this is the case, the method of increments is ac¬ 
tually retrieving more correlation energy than DMRG (as 
one can see the error is negative in a quite wide regime). 
Only in the metal-like regime this approach gives worse 
results with respect to the other reported even if it un¬ 
derestimate the module of E corl only by 2 x 10~ 3 Eh with 
respect to RAS(4,24)). This can be explained considering 
the high delocalization of the wave function in this regime 
which is hard to describe with a local approach. A quan¬ 
titative comparison between the results with CAS-Mol, 
DMRG, CCSD(T) and different CAS and RAS methods 
is reported in Table [TJ Again one can realize the difficulty 
in retrieving a meaningful amount of the correlation en¬ 
ergy with canonical methods and how in general the local 
approach gives reasonable results with relatively less ef¬ 
fort. 

Let us now consider the strong deviations shown in Fig. [3] 
(dashed lines) that occur when the CAS-Mol is employed 
starting from a HF configuration which is not the dom¬ 
inant one. In order to explain why in this situation the 
method fails in describing the electronic structure of the 
system, let us analyze the behavior of the individual in¬ 
crements as shown in Fig. [Gj What can be observed is 
that, as we overcome the crossing, the required conver¬ 
gence stated in Eq. [5] and [6] break down and higher or¬ 
der increments give comparable contributions to E co „ as 
lower order ones. This is equivalent to say that in the em¬ 
ployed orbital basis (i.e HF configuration), higher order 
contributions become more and more important which 
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is analogous to what was concluded from the analysis 
of the mutual information as described in our DMRG 
analysis^ that showed the increase of long range entan¬ 
glement. This strong dependence on the starting configu¬ 
ration is of course a great disadvantage of the Mol which 
can otherwise yield accurate results with relatively cheap 
calculations. 

Since we are forced to calculate both configurations, it 
might appear that no particular advantage arises from 
the use the CAS-Mol with respect to the single reference 
approach. On account of this, in Fig. [7] we report the 
differences between the energies obtained from Conf 1 
and Conf 2 in the crossing regime for both methods. As 
it can be seen, in the case of CAS-Mol this difference is 
around 10 mE/, smaller than CCSD(T) and ranges be¬ 
tween ±5 m Eh. This means that in this narrow regime, 
the two bases give us comparable results leading to a 
smooth curve, which is not the case for CCSD(T). In ac¬ 
count of this and previous statements, we can conclude 
that the CAS-Mol allows to describe the behavior of the 
PES also in the crossing regime allowing us to describe 
the whole dissociation curve. 

We conclude this section comparing the individual in¬ 
crements as obtained from CCSD(T) and CAS-SCF (see 
Fig. [G] and Table HI]) . For the former method, the lower 
order increments are larger than for CAS-Mol because a 
larger active space is involved, but at the 3- and 4-body 
level the orbital optimization which retrieves the static 
correlation starts playing a bigger role than the size of 
active space. In general the single reference method con¬ 
verges faster, but the total correlation energy is larger for 
CAS-SCF than for the single reference approach. 



Be-Be (A) 

FIG. 4. (Color online) Potential energy curve for Be6 calcu¬ 
lated using the Mol at the 4-body level within the CAS-SCF 
approach starting from the two HF main configurations. 


B. Extension to larger rings 

After analyzing the use of the method for Beg we report 
the results for larger rings in order to evaluate the conver- 
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FIG. 5. (Color online) Comparison between the differences 
of the PES for Beg obtained through different methods and 
DMRG(M = 1024). 


gence towards the thermodynamic limit. DMRG allowed 
us to describe the dissociation of Beio (20 electrons in 40 
active orbitals) and as for the smaller cluster we report 
(Fig. [8J) the energy differences of CAS-Mol and CCSD(T) 
with respect to DMRG. Also in this case it can be seen 
that the former method behaves much better than the 
single reference one and that the obtained values differ 
from DMRG ones by around 5 m Eh- 
As already stated, moving forward the aim of describing 
periodic systems, local methods become better candidate 
than canonical ones for obtaining accurate correlation 
energies. Within the CAS-Mol, the amount of necessary 
CSFs depends on the body order only and not on the size 
of the system, while in a canonical approach prohibitive 
active space sizes become necessary. The PES calculated 
for Beio and Bei 4 rings at the 4-body level are shown in 
Fig. |9] in comparison with Beg. Even from these small 
cluster we can deduce the trend for increasing size, that 
is a faster convergence in the insulator-like regime than 
in the metal-like regime and a shift of the crossing to¬ 
wards larger internuclear distances. 

In order to highlight these convergences we report in Ta¬ 
ble eh the individual increments and the calculated cor¬ 
relation energies for Be n rings with n = 6, 10, 14, 22, 30 
and 90 in the two regimes. As one can see in Fig. fTOlthe 
individual increments can be fitted with a n~ u function. 
This way through extrapolation it was possible to eval¬ 
uate the values for n —> oo, i.e. for the periodic chain, 
which are also reported in Table Hill 
It can be observed that both in the metal-like and 
insulator-like regime the trend is basically quadratic for 
the 1-body while it rapidally deviates for higher order in¬ 
crements in the metal-like regime where the depencende 
on the angle is stronger. It has to be underlined that for 
the 4-body case, the value for n = 6 was not included 
in the fit because for obvious geometrical reason was to 
much out of trend. Clearly the 1-body dominates the 
total correlation energy and therefore E COII depends also 
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FIG. 6. (Color online) Percent of the correlation energy contribution of every incremental order as a function of distance for a 
Be6 ring. The data are reported for configuration 1 (doubly occupied sp hybrids) and configuration 2 (2s 2 ). The convergence 
criteria fulfill only in particular distance regimes where the incremental scheme can successfully be used and more than 99% of 
the correlation energy can be retrieved at the 4-body level. At the crossing point 2.60A, both configurations can be used. 
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Fock energy too. 

Of course a faster convergence would be achieved if lin¬ 
ear chains had been considered since no angle depen¬ 
dence would be present as in the ring shaped cluster, 
but these were chosen to impose a higher symmetry to 
the system, deleting any border effect and allowing the 
different bodies to be equivalent to each other. More¬ 
over the cyclic structure recalls more the periodic Born- 
von-Karman boundary conditions imposed in the peri¬ 
odic case. 


V. CONCLUSION 


FIG. 7. (Color online) Comparison between the energy dif¬ 
ference F/Confi — Econi 2 for the canonical CCSD(T) and the 
CAS-Mol around the crossing for Be6- At the Full-CI level 
this difference is expected to be zero. At the 4-body level the 
method of increments gives an acceptable result. 


quadratically on 1/n. We conclude showing in Fig. [Tl] 
that a similar behavior was obtained for the Hartree- 


The method of increments has been applied using a 
multireference formalism in order to describe the disso¬ 
ciation of ring shaped beryllium cluster of different sizes. 
The large static correlation involved in these systems 
made them good candidates for our analysis. Being size 
consistent and size extensive the Mol allowed us to ob¬ 
tain CAS-SCF wave function that can be further used 
as basis for multireference calculations. Key point of this 
investigation was that we were able to describe the whole 
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TABLE I. (Color online) Comparison between different canonical CAS and RAS approaches, CCSD(T), DMRG calculations 
reported in Ref.[l9(| and Mol calculation at the 4-body level as described in the present paper for Be6. 



Metal-like - 2.10A- 

Conf 1 

Insulator-like - 3.00A- 

Conf 2 

Calculation 

CSEs Energy (hJ h ) F/’ C orr(%) 

CSEs Energy ( E h ) 

A'corr(%) 

CAS(6,6) 

34 -87.329677 

8.09% 

34 -86.947385 

8.54% 

CAS(6,12) 

2,086 -87.364902 

29.31% 

2,086 -87.034032 

33.57% 

CAS(6,15) 

8,155 -87.367285 

30.74% 

8,155 -87.035797 

34.08% 

RAS(4,21) 

13,523 -87.373196 

34.30% 

13,521 -87.002608 

24.49% 

CAS(6,21) 

64,835 -87.373263 

34.34% 

64,835 -87.002788 

24.54% 

RAS(6,18) 

1,268,308 -87.425535 

65.82% 

1,268,308 -87.146271 

65.98% 

RAS(4,23) 

134,843 -87.428804 

67.79% 

134,837 -87.141082 

64.48% 

RAS(2,24) 

884 -87.455551 

83.90% 

885 -87.160246 

70.02% 

CCSD(T) 

-87.479226 

98.15% 

-87.261723 

99.32% 

RAS(4,24) 

295,746 -87.479826 

98.51% 

295,752 -87.240226 

93.12% 

RAS(6,24) 

14,972,954 -87.482727 

100.26% 

14,972,954 -87.260875 

99.08% 

CCSD(T)-MoI(4-body) 

-87.479272 

98.18% 

-87.261705 

99.32% 

MR-MoI(4-body) 

866,320 -87.478599 

97.77% 

866,320 -87.264183 

100.03% 

DMRG(M = 1024) 

-87.482294 

100.00% 

-87.264065 

100.00% 


TABLE II. Individual incremental orders and total correlation energy of Be6 calculated with the Mol at the 4-body level using 
a single reference method, CCSD(T) and the multireference method as described in this paper. 


Configuration 1 Configuration 2 



Metal-like - 2.10A 

Crossing regime - 2.60A 

Crossing regime - 2.60A 

Insulator-like - 3.00A 

Increments 

CCSD(T) CAS-SCF 

CCSD(T) 

CAS-SCF 

CCSD(T) 

CAS-SCF 

CCSD(T) CAS-SCF 

Ei 

Ei,,- A Uj 

Ei.i.fc.i A tijkl 

-0.0197077 -0.0195083 

-0.0072719 -0.0051364 

-0.0001276 -0.0019666 

-0.0000652 -0.0004490 

-0.0282556 

-0.0109604 

-0.0004272 

-0.0000009 

-0.0278839 

-0.0088279 

-0.0043662 

-0.0007402 

-0.0443223 

-0.0113016 

-0.0005745 

-0.0001604 

-0.0438349 

-0.0105205 

-0.0028508 

-0.0001981 

-0.0532219 -0.0531172 

-0.0040539 -0.0041013 

-0.0000447 -0.0004985 

0.0000035 -0.0000157 

F, 

J-'corr 

-0.0271724 -0.0270603 

-0.0396440 

-0.0418181 

-0.0563588 

-0.0574042 

-0.0573169 -0.0577327 

Etot 

-87.479272 -87.478599 

-87.300283 

-87.3133272 

-87.310961 

-87.317234 

-87.261705 -87.264200 


TABLE III. Individual incremental orders calculated with the CAS-Mol for Be„ ring-shaped cluster of different sizes form n ■ 
6 to 90 for two internuclear distances, 2.10A and 3.00A. The extrapolated values for the infinite chains are also reported. 

Metal-like regime - 2.10A Conf 1 

Increments 6 10 14 22 30 90 n —¥ oo 

Ei 


Ei,i Ae»j 


-0.01950832 -0.01852919 -0.01826832 -0.01810244 -0.01804935 -0.01799396 -0.017990(2) 
h] —.j -0.00513636 -0.00388669 -0.00358258 -0.00340908 -0.00335721 -0.00330521 -0.003305(3) 

Ei.j.fc A ^ijk -0.00196657 -0.00152785 -0.00145735 -0.00141879 -0.00140766 -0.00139659 -0.001402(4) 

]T. . ; A e ijkl -0.00044904 -0.00048085 -0.00047395 -0.00047172 -0.00047140 -0.00047112 -0.00047118(4) 

£corr -0.02706030 -0.02442458 -0.02378220 -0.02340203 -0.02328562 -0.02316688 -0.023155(2) 


Increments 


Insulator-like regime - 3.00A Conf 2 
10 14 22 30 


90 


Ei U -0.05311719 -0.05218930 -0.05190679 -0.05174649 -0.05169682 -0.05164440 -0.051633(7) 

Eij A tij -0.00410130 -0.00450650 -0.00464482 -0.00472404 -0.00474799 -0.00477369 -0.004783(6) 

Ei,,- k A ^ijk -0.00049846 -0.00063636 -0.00068277 -0.00070780 -0.00071598 -0.00072474 -0.000727(2) 

Ei,,-, fc ,i Aejjfci -0.00001574 0.00002165 0.00002305 0.00002433 0.00002457 0.00002476 0.0000247(5) 

Bcorr -0.05773268 -0.05731051 -0.05721133 -0.05715400 -0.05713622 -0.05711807 -0.057119(2) 




















Be-Be (A) 


FIG. 8. (Color online) Energy difference between CAS-Mol 
at the 4-body level and DMRG(M = 1024) for the Beio- 
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FIG. 9. (Color online) Potential energy surfaces for Be6, Beio 
and Bei 4 rings obtained through the CAS-Mol at the 4-body 
level. 


PES also where single reference methods fail, i. e. close to 
the crossing. DMRG calculations previously performed 
were crucial to evaluate the reliability of our results. The 
correlation energy of a system as large as Bego was ob¬ 
tained using this approach. We want to underline that 
within a canonical method this calculation would involve 
180 active electrons in 360 active orbitals which without 
a local approach this would be impossible to treat. Fi¬ 
nally, investigating different sizes of Be„ rings we could 
evaluate the correlation energy for the periodic system. 
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Metal-like regime - 2.10A Conf 1 
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Insulator-like regime - 3.00A Conf 2 
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FIG. 10. (Color online) Convergence with the ring size of the largest individual increments of each order. The red lines were 
obtained by fitting the calculated data with an equation of the form a/n v + E°° where E°° is the value extrapolated for the 
periodic chain. 
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FIG. 11. (Color online) Convergence with the ring size of 
3.00A(right). 
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Hartree-Foc.k and total energy per atom for 2.10A(left) and 
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